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ABSTRACT 

We study gravitational lensing statistics, matter power spectra and the angular power 
spectra of the cosmic microwave background (CMB) radiation in x-matter models. 
We adopt an equation-of-state of x-matter which can express a wide range of matter 
from the pressureless dust to the cosmological constant. A new ingredient in this 
model is the sound speed of the x-component in addition to the equation-of-state 
Wo — Pxo/PxO- Except for the cosmological constant case, the perturbations of x- 
matter itself are considered. Our primary interest is in the effect of non-zero sound 
speed on the structure formation and the CMB spectra. It is found that there exist 
parameter ranges where x-matter models are consistent with all current observations. 
The x-matter generally leaves imprints in the CMB anisotropy and the matter power 
spectrum, which should be detectable in future observations. 

Key words: cosmology: theory - dark matter - gravitational lensing ~ large-scale 
structure of the universe - cosmic microwave background 



1 INTRODUCTION 

There is growing observational evidence that the mean mass density of the universe is less than the critical density. The growth 
of large scale structure in cold dark matter(CDM) models requires F — Q,Mh = 0.25±0.05 (Peacock & Dodds 1994) where Q.m 
is the ratio of the matter energy density to the critical density and h is the Hubble parameter in units of lOOkms^^Mpc" 



Recent observations of high-redshift Type la supernovae also indicate 51a/ < 1 (Garnavich et al. 1998 



199^ 



Pcrlmutter et al 



Schmidt et al. 1998; Riess et al. 1995). The observations of the age of the globular clusters 12 ± 2Gyr (Chaboyer et al 



1995) and the Hubble parameter Hq = 65 ± lOkm/s/Mpc (Freedman 1996; Riess et al. 1996), though there might be some 



modifications by recent Hipparcos results (Feast & Catchpole 1997; R,eid 1997), indicate HqIq ^ 0.8 where to is the age of the 



universe. All these observational pieces of evidence may not be explained in the framework of standard CDM model(sCDM) 
in which Q.m ~ 1- The non-zero cosmological constant (ACDM)may be needed to increase the age of the universe as well 



as to make the total energy density equal to the critical density as generically predicted by inflationary cosmology (Ostriker 



& Steinhardt 1995). However, as noted in Nakamura et al.(1995), at least for the age problem, the cosmological constant is 



not the only way to solve the problem: if Hoto > 1, then what is needed in general is the "matter" which violates the strong 
energy condition such that p + 3p < 0. 

Turner and White (1997) proposed xCDM as an alternative to ACDM models in which a nonclustering matter of px ~ wp^ 
{w being constant with — 1 < ui < 0) resides in addition to ordinary CDM so that the universe has a flat geometry. They 
also noted some topological defects behave like the "x-component". However, as they admit, a component with < is 
highly unstable to growth of perturbations on small scales. Therefore they assumed xCDM to be a smooth component on 
small scales. The notion of a smooth component is ill-defined because i) the smoothness is not gauge-invariant and ii) it is 
unphysica l to ignore the respo nse of the component to the inhomogeneities in other cosmological fluid under the equivalence 
principle (Caldwell et al. 1998). 



Recently, we have constructed a phenomenological energy momentum tensor for "x-matter" ( []hiba et al. 1997 ) which i) 
is locally stable, and ii) obeys the causality. To satisfy the local stability requirement, we have to separate the sound speed 
Cj = p'x/px from the coefficient of the equation of state = p^j p^. We assume x-matter obeys the equation of state such that 
Px ~ Wo pxo 4- (px — PxO ) , where px , Px and pxo are the pressure, the energy density and the present energy density of x-matter. 
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respectively. As regards the background universe, this x-matter model contains both Einstein-de Sitter {wo = — 0) and 
A universes {wo — —1) as a limit. A new ingredient in this model is the sound speed of the x-component in addition to the 
equation-of-state wo- Since at present we do not know how the vacuum energy interacts gravitationally (Weinberg 1989), a 
phenomenological approach taken here is worthwhile. 

On the other hand, when we deal with perturbations of x-matter to study the structure formation in the universe, the 
limit of the A universe in which density and pressure stay constant has to be carefully taken. In this case, constant density 
and pressure allow the existence of only smoothed distribution from the gauge-invariance. Therefore x-matter is equivalent to 
the cosmological term. 

In this paper, we shall make further observational tests of x-matter models as alternatives to sCDM or ACDM. First we 
consider the statistics of gravitational lensing of quasars. Then we study the evolution of density perturbations of x-matter 
models. Accordingly, we obtain observational quantities such as the matter power spectrum and cosmic microwave background 
(CMB) anisotropies. The class of cosmological models we consider are spatially flat Friedmann-Robertson-Walker(FRW) 
models which are generically predicted by cosmological inflation. We consider x-matter, dust matter (baryon and/or CDM) 
and radiation (photon and neutrino) while we do not include the cosmological constant term whose density is smoothly 
distributed. Similar analysis was made in xCDM models (Bloomfield Torres & Waga 199f ; Silveira & Waga 1997), in scaling 



field model (Ferreira & Joyce 1997a; Ferreira & Joyce 1997b), the decaying lambda model (Viana & Liddle 1998) and in 
quintessence(QCDM) model ( Caldwell et al. 1998| ). When our paper was at the nearly end of completion, we became aware 
of the recent paper by Hu (1998) where similar models were considered. In his paper, he employed very general equation of 
states which include our model. However, since we focus on the simple, yet rather general equation of state, we could perform 
thorough investigation and observational tests of the model here. 

The lensing frequencies are expected to be reduced in x-matter models, since in a x-matter model the distance to an 
object at redshift z is smaller than the distance to the same object in A-dominated models. We study the lensing frequencies 
and set quantitative limits on the present energy density of x-matter using quasar surveys. 

The height and location of the peak in the angular power spectrum of CMB depend on the matter-radiation equality 
and the expansion rate of the universe at the surface of last scatter. Different equation-of-states yields the different matter- 
radiation equality time, which would leave imprint in the peak height. Non-zero Cs affects the expansion rate of the universe, 
and its effect would leave the imprint in the location of the peak even if for low fix. 

This paper is organized as follows. In §2, we review our x-matter model. In §3, we present calculations of gravitational 
lensing in this model. In §4, a detailed perturbation theory in x-matter models is given and the spectra of cosmic microwave 
background are calculated. §5 is devoted to summary. We use the units of c = 1. 



2 X-MATTER MODEL 



In this section we review our x-matter model (Chiba et al. 1997) 



x-matter possesses following nature, i.e., Wx = Px/px < to elongate the cosmic age and = Sp/Sp > to stabilize 
growth of perturbations. Cs represents the sound speed of x-matter and determines how the density evolves. Together with 
the latter condition, we demand that it is not faster than the speed of light from the causality. As a simple form, we employ 
the equation of state of x-matter as 

Px = IfOPxO + Cs(px - Pxo), (1) 

where the suffix denotes the present value. Although wo and Cj, are functions of px in general, to simplify the discussion we 
assume that they are constants. Our x-matter model includes the cosmological constant model as a special case since px ~ — PxO 
for Wo — —1 regardless of the value of cj. However, as is mentioned in the previous section, it should be noticed that x-matter 
is equivalent to the cosmological term and it is no longer 'matter' in this limit. According to above requirements, we consider 
the region —1 < wo < 0, < < 1, and px > 0. From the requirement that the process of Big Bang Nucleosynthesis (BBN) 



is not appreciably disturbed by the existence of x-matter, is further constrained: < 0.15 (Chiba et al. 1997). 
The energy equation is 

where a is the scale factor. Normalising a as ao = 1, we have 

PxO /(-i I \ -3(l+c^) , 2 \ /Q\ 

l^d V ' +CB-woj. (3) 

Note that uio = — 1 corresponds to the cosmological constant model irrespective of c^. We also note that the cosmological 
constant model is an asymptotic limit in our x-matter model: px/px —1 as a ^ oo. The effective scalar field potential for 
the x-matter model is described in Appendix A. 
The Friedmann equation in a flat universe is 



Px 
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where ^m, fir and Q,^ axe the present density parameter of matter(baryon and CDM; Q.m — fie + f^CDM), radiation (photon 
and neutrino) and x-matter, respectively and over dot denotes the time derivative. 
The cosmic age is given by 

Hoto ^ J^'^ (ilMa'-' + ^ ((1 + u,o)a-^(^+==> + - wo) ) (5) 

Here we have neglected the energy density of the radiation components since the effect of them on the cosmic age is negligible. 
The functional form of Hgto shows that it is a decreasing function of Cs or wq. In Fig.l, the contour plots for Hoto in the 
{wo, fix) plane are shown {flm = 1 — i^x). 



3 GRAVITATIONAL LENSING 



Now we consider the statistics of gravitational lensing. We compute the predicted number of gravitational lensing events 
taking into account the magnification bias and compare it with the observational data. Our analysis largely depends on 
Kochanek (1996) since our main interest is in the comparison with A-dominated models. Accordingly, we will neglect the dust 
obscuration effect considered by Malhotta et al.( 1997 ) and the evolution effect of lensing galaxies which would weaken the 
constraint. 



3.1 Lensing Probability 



We consider a singular isothermal model for a lens mass distribution mainly for simplicity. Hubble Space Telescope (H ST) 
observations show that E/SO galaxies are effectively singular (Tremaine et al. 1994), which may support our assumption. The 
effect of a finite core radius is shown to be not so significant (Kochanek 1996) because a core radius increases the cross section 



while decreases the amplification compared to a singular isothermal model and both effects are compensated in the lensing 



probability. The cross section for lensing events for a singular isothermal sphere lens model is given by (Turner et al. 1984) 
(see Appendix B for its derivation) 



1^34/ DolDls 
167r V — J- 

V Dos 



(6) 



where i; is the velocity dispersion of the dark halo of a lensing galaxy, Dol is the angular diameter distance to the lens, Dos 
is the angular diameter distance to the source, and Dls is the angular diameter distance between the lens and the source. 
The differential probability of a lensing event is 

dt 



dr 



no{l + zl) Ba^f—dzL, 

dZL 



dt 



= IGtt no« (1 + 2:l) B\ — - — dz-L, 

\ Dos ) dzL 



(7) 



where no is the present-day number density of lensing galaxies, and B is the amplification bias factor defined below. The total 
lensing probability is obtained by integrating Eq.(^ with respect to zl- 



The number count of galaxies is fitted by a Schechter function (Schchter 1976) 

-L/L,dL 

(?!>. = 0.026 ± 0.08/l^Mpc-^ a 



4,g{L)dL = 4>, 



(8) 



Recent Autofib survey gives 

from E/SO type galaxies. The contribution from spiral galaxies is not so important (Kochanek 1996; Fukugita & Turner 



1.09lo og (Ellis et al. 199e|). We only consider the contribution 



1991 



Maoz & Rix 1993) because the typical velocity dispersion (c 

3 



140km/s) is small. The comoving number density of E/SO 



"=*^^Mpc- 



(Heyl et al. 1997) 



galaxies Ue is rie — 3.5 x 10 

We assume that the velocity dispersion of the dark halo v is related to the luminosity by a "Faber- Jackson relation" of the 
form 11 = 11, (L/L,)^/^. For a singular isothermal model, Kochanek found that v, = 225 ± lOkm/s and 7 = 4.1 ± 0.9(Kochanek 
1994). Following Kochanek (1996), we do not consider the (3/2)^''^ correction in the velocity dispersion since the dynamically 



estimated velocity dispersion of dark matter is not different from the central velocity dispersion. 

As for the distance formula, we adopt the standard angular diameter distance dA in a FRW universe given by 

'1/2 



dAHo 



1 



du 



l + ci 



(' 



l + wo){l + u) 



3(i+c;) 



- Wo 



(9) 



Then the angular diameter distance between the lens and the source Dls is given by 
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DlsHo = 



l+ZS 



du 



nM(l + uf + ((1 + Wo){l + uf ^+'^i^ + - ™o) 



-1/2 



(10) 



The magnification bias B is an enhancement of probability that a quasar is lensed. The bias for a quasar at redshift z 
with apparent magnitude mg is written as 

/ druQ (jiQiniQ +2.5log A, Zg)P{A)dA 



B{mQ,ZQ) = 



(11) 



/ dmQ<j>Q{mQ,ZQ) 

where <j)Q{mQ, zq) is the luminosity function of quasars at redshift zq, Amin is the minimum total flux amplification, and 
P{A) is the probability distribution for a greater amplification A. 

We use the Kochanek's "best model" ( Kochanek 1996) as a quasar luminosity function. In this model, like Boyle, Shanks 
and Peterson's model ( Boyle et al. 198^ ), the luminosity function has the form 



(12) 



where the bright -end slope qq and faint-end slope /3q are constants, while the break magnitude mo{z) evolves with redshift 
( iKochanek 1996| ) 

mo + {z - 1.0), z <l 

mo{z) — ^ mo, 1 < 2 < 3 (13) 

mo — 0.7(2 — 3), z > 3. 



Fitting this evolution model to the quasar luminosity function data in Hartwich and Schade (199C) for z > I, Kochanek find 
that "the best model" has aq = 1.07 ± 0.07, /3q = 0.27 ± 0.07 and mo = 18.92 ± 0.1 6 at B magnitudes. He noted that this 
model is no better than Wallington and Narayan model (Wallington & Narayan 1993). 

For a singular isothermal sphere, the amplitude amplification A{r) from a projected lens-source separation r is 



A = 2bcr/r, 



(14) 



where bcr is the critical impact parameter given in Eq.(B6). Then A> 2 for a singular isothermal sphere. The probability of 
amplification is hence 

l3 



P(A) = 8/A^ 

Putting all these pieces together, the lensing probability is finally given by 

n, dt .3 r dL if LY -L/L,,^ :i ( L\^/^ , ( DqlDls 



(15) 



dzr 



dt 



' H-^dzL 



(1 + ZLfB{mQ,ZQ) 



( DolDl 



Here F is given by 



F = 16n^ Ho ■'n,vtT{4:/-/ + a + 1) = l.GZo s x 10 



(16) 



(17) 



where r(a;) is the gamma function. F roughly means the probability that a quasar at horizon distance is lensed by a L* 
galaxy. 

Fig. 2 shows the optical depth in a flat universe as a function of a source redshift Zs for Q]\f — 0.1, 0.3, 0.5, 0.7, respectively. 
Current constraint (Kochanek 1996) may be translated as t{z3)/t{zs;Qm = 1) ^ 5. We expect that x-matter dominated 
models may survive the observational test. Next we estimate the allowed parameter regions more quantitatively. 



3.2 Comparison with Observations 



We fix the parameters in the quasar luminosity function (ag , I3q , mo) but allow changes of parameters in the galaxy luminosity 
function within their error bars. We then compare the predicted expectation number of lenses with the observed one. Data 



are taken from HST snapshot survey (Maoz et al. 1993) and the NOT survey (Jaunsen et al. 1995). We use only quasars at 
z > 1. The total number of quasars in the sample is then 625. The sample contains 4 lensing events caused by galaxies(0142- 
100,1115-1-080,1413+117,1268+1011). Since most of the quasars in the sample do not have a B-band magnitude, we transform 
mv to a B-band magnitude using B — V — 0.2 as suggested by Bahcall et al.(Bahcall et al. 1992). Fig. 3 shows the expected 
number of lenses for wo — —1, —0.8, ^0.6, —0.4, ^0.2, and = 0,0.05,0.1,0.15 as a function of Jlx. Each curve from 
top to bottom corresponds to wo = —1, —0.8, —0.6, —0.4, —0.2,0. wo = —1 corresponds to the cosmological constant model 
irrespective of cj. The expected number increases as fix increases. However, its rate decreases for increasing wo or c^. 

We perform a maximum-likelihood analysis to determine the confidence level of model parameters. The likelihood function 

is 
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(18) 



where Nu is the number of unlensed quasars, Nl is the number of lensed quasars, and Pi is the probabihty that a quasar 
is lensed. L is a function of three parameters (c^,wo,^^x)- We allow parameters and wq to vary in the range < < 
0.15, —1 < Wo < 0, respectively. The logarithm of the ratio of the likelihood to its maximum —21x1 L / Lmax is asymptotically 
distributed like a distribution with the degrees of freedom of parameters involved (Kendall & Stuart 1973). Thus we can 
determine the confidence level. 

Likelihood contours are plotted in Fig. 4 projected onto constant Cs plane plane. Likelihood levels are, from left top to 
right down, 68.3%, 95%, 99%. We find that larger f2x is allowed as increasing uiq. Changing Cs has little effect on likelihood 
contours because the effect of can be significant only for much higher redshift. The 95%C.L. level can be well fitted 
by f2x < 0.85 + (1 — 4Cs/3)(u)o + 1). The constraint on the cosmological constant model(iuo ~ —1) is less stringent than 
Kochanek's because of the difference in the estimated number density of E/SO type galaxies (Ellis et al. 1996). Taking the 
same number density of galaxies as Kochanek's, we obtained the constraint Ao < Q.67(95%C.L.) for flat models consistent 
with his constraint; Aq < 0.66 (95%C.L.). 



4 PERTURBATIONS IN X-MATTER MODELS 
4.1 Evolution and General Features 



The gauge invariant perturbation equations for generalized fluid, i.e., x-matter are given in the Fourier space as (Kodama & 



Sasaki 1984) 



(T^)' = -fcK + -^^iVcip.A. , (19) 

Vl+Wx/ l+WxCSp ^ 

a 

and 

K' = --(1 - 3c?)V; + f-^Ax - ^-v] + fc* , (20) 
a \ 1 + Wx a y 

where Ax, axe density and velocity perturbations of x-matter, k is the wave number, and = p^/p^ are sound speed 
and pressure-density ratio of x-matter, ' is the derivative in terms of conformal time rj, c^, pa and Ac are sound speed, 
density, density perturbations of the a component, respectively, indicates summation in terms of all components, V is 
total velocity perturbations and 'i' is the gravitational potential. Above density perturbations are ones relative to the total 
matter rest frame. We solve these equations numerically together with equations for ordinary baryonic matter, radiation 
(photons and neutrinos) and/or cold dark matter. The time evolutions of x-matter density fluctuations are shown in Figure 
5. In this figure, we take h = 0.5, Qb ~ 0.05, Six = 0.95 and Q,cdm ~ 0. For the comparison, we plotted the evolution of CDM 
perturbations (a dotted line). 

If we set Cs to be zero in order to study the dependence on wo, we find decay of fiuctuations near the present epoch (a 
dashed line). It is interesting to examine what causes this turnover. In case of CDM dominated models, we can assume 'i' is 
constant in time, and the scale factor a ccrj^ . Therefore from equations ( p^ ) and (po|), we obtain ocq and Ax oc 77^(1 + Wx)- 
Accordingly, we find Ax ex ?7(1 -I- Wx)(3Wx-|- 1). Thus, we can conclude that Ax decays if Wx < —1/3 in case of CDM dominated 
models. This critical value Wx = ~l/3 corresponds to the breaking of the strong energy condition. For x-matter dominated 
models, we numerically find a slightly larger critical value (wx — —0.28). It is because even for higher in these models, 
decay of the self-gravitational potential makes perturbations difficult to grow. 

When we take a non-zero value of c^, on the other hand, we find the acoustic oscillation of fluctuations of x-matter after 
crossing the Jeans scale. In Figure 6, the Jeans mass scales for various f2x and are shown as a function of redshift. The 
Jeans mass of x-matter is defined as 
An / A - ^ 3 



47r f X-iY , N 



where 

/ \ 1/2 

Here ptot is the total density defined as ptot = Px + pcDM + ps. Solid and dashed lines denote wo = —0.2 and —0.8, respectively. 
In this figure, we find that it is difficult to form the structure below galaxy scales if > 10~* when x-matter is the dominant 
component of the universe. However, the constraint can be loosened if we allow the time variation of c^. 

To see the effect of wq and Jeans oscillation, the matter transfer functions are shown in Figure 7. In panel (a), we take 
Cj = 0. In this panel, smaller wq indicates smaller amount of x-matter in past and later matter-radiation equality since we 
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fix the present x-matter density. Therefore the location of the knee in the transfer function which corresponds to the horizon 
scale of the matter-radiation equality shifts to the larger scale when we take smaller wq. In panel (b), on the other hand, we 
fix Wo to be and change c^. The Jeans oscillation is seen in the transfer function even in case of ci = 10"^ 

In Figure 8, the angular power spectra of CMB anisotropies, so-called Ce are shown for the same models as Figure 5. The 
shape of Ce for the — model is quite similar to the one with the cosmological constant. The peak locations and the shape 
on small £ (large scales) of this model are different from sCDM which is the fiducial CDM models with h = 0.5, f2o = 1 and 
Qb ~ 0.05. A dent on the large scale {I ~ 15) is caused by the decay of gravitational potential just like lambda models (see 
e.g., Sugiyama & Silk 1994). Since the volume of the universe of the negative wq model is bigger than the one of the standard 
Einstein-de Sitter model, the peak locations shift to the smaller scales as is shown in this figure. 

On the other hand, Ce of the model with high has a quite different shape. In this model, the universe has never become 
purely matter dominated or radiation dominated. In other words, the gravitational potential has been decaying all the time. 
This decay boosts Ce in the long range of ^'s. 

One interesting feature of this figure is that the damping locations at high £ of all Ce's are almost the same. It is simply 
because recombination history is the same regardless of the equation of state since it is determined by the ratio between the 
number densities of photons and baryons. 



4.2 Observational Constraints 

Now let us make a comparison with observational data. In Figures 9 (/i = 0.5) and 10 (h — 0.7), we plot the contour maps 
of (78, which is the rms density field in a sphere of radius 8/i~^Mpc with the top hat window function, on the wq — fix 



plane. Here we employ COBE 4yr normalization (Bennett et al. 1996; Bunn & White 1997). The observational value of as 
is as = (0.52 ± 0.04:)Qq'^^^'^'^^^° for fiat models from the analysis of the local cluster X-ray temperature function (Eke et 



al. 1996 ) . Here flp is the density parameter of total matter in case of fiat A-models. Similar results are obtained from other 



analyses ( White et al. 199S ; Viana fc Liddle 1996 ). While there is Qo dependence in the observational value of as, we expect 
the dependence of as on fix and wo is relatively weak unless wq — —1 and fix — 1. 

Employing erg ~ 0.6, we plotted total matter transfer functions in Figures 11 (panels (a)-(c) for h — 0.5 and panels 
(d)-(f) for h — 0.7) for various Cs's and wo's. In case of < 10~^, transfer functions do not depend on since the Jeans 
scale is smaller than lh~^Mpc even at present. Therefore the shape of transfer functions is merely controlled by the matter- 
radiation equality epoch since the horizon scale at this epoch corresponds to the knee of the transfer function. Regarding 
the normalization of perturbations, CMB anisotropies are affected by the decay of gravitational potential, which is so-called 
integrated Sachs- Wolfe effect. Hence one might expect different normalization factors from the same shape of the transfer 
function for models with different wq's. However, the decay of the gravitational potential caused by wq takes place at very late 
epoch even in case wo — —1 as is shown in Figure 5. Therefore the integrated Sachs- Wolfe effect does not provide dominant 
contribution on CMB anisotropies on the COBE scale unless the universe is dominated by x-matter. Eventually we conclude 
that fixing as indicates fixing the matter-radiation equality and fixing the shape of the transfer function. On the other hand, 
we find a mild damping feature in the transfer function on small scales just like a mixed dark matter model(i.e. cold+hot 
dark matter) if we take > 10"** because of the Jeans oscillation of x-matter perturbations. In particular, it is interesting 
that models with fix — 0.1, — 0.15 and wo = —0.8 fit fairly well to the observational result by Peacock and Dodds (1994) 



as shown in Figure 12. 

In Figure 13 {h = 0.5 for panels (a)-(c) and h = 0.7 for panels (d)-(f)), we plot Ce's of models with as — 0.6. In case of 
Cg < 10~®, Ce's do not depend on neither Cs, wq nor fix (see panel (a) or (d)). Comparing to the pure CDM model, however, 
we find higher acoustic peaks. It is because the epoch of the matter-radiation equality is the same for models with the same 
erg. On the other hand, the equality epochs of the pure CDM models whose erg are 1.2 and 1.6 for h = 0.5 and 0.7, respectively. 



are earlier than ones of x-matter models. It is well known that the later the equality epoch, the higher the peaks are (Hu et 



al. 1997). It is the same for panels (b) or (e). Smaller wo indicates smaller amount of x-matter in past and the later matter 
radiation equality since we fix the present x-matter density. Therefore it is shown higher peaks for smaller wq. In case of larger 
Cs, on the other hand, the peak locations are shifted on smaller scales (larger £) as is shown in panels (c) or (f). With larger 
Cs, the equation-of-state of x-matter is in between radiation and matter and x-matter is a (sub)-dominant component in the 
universe around the matter-radiation equality epoch even for small fix. This component causes the decay of the gravitational 
potential since only pure radiation or matter domination keeps the potential staying constant. And this decay boosts the first 
acoustic peak as is shown in panels (c) or (f). 

In conclusion, adopting erg ~ 0.6 constrains x-rnattcr content fix ^ 0.1 for ~ 0.1, or if we consider the universe to be 
x-matter dominated, then its sound speed is constrained as < 10~® at least at the epoch of galaxy formation. Former case 
can fit the observational data by Peacock-Dodds well. With high precision measurements expected from MAP0(Microwave 



* http://map.gsfc.nana.gov 
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Anisotropy Probe), PLANCkQ, 2Df|(2 Degree Field) and SDSs|(Sloan Digital Sky Survey), we should distinguish x-matter 
models from sCDM or ACDM. 



5 SUMMARY 

We have investigated the possibility that the dark matter component has the equation of state px = wopxo + <?s{px — pxo) 
such that — 1 < uiQ < — 1/3 and < 0.15. We have studied the detailed observational consequences of x-matter models 
including the statistics of gravitational lensing of quasars, matter power spectrum and cosmic microwave background (CMB) 
anisotropies together with a comprehensive study of density perturbations of x-matter models. The gravitational lensing 
statistics constrains Q,^ as Q,^ < 0.85 + (1 — 4:Cg/3){l + it)o)(95%C.L.). We have shown that introducing the sound speed of 
x-component has a great effect on CMB anisotropies and matter power spectrum. We find that adopting erg ~ 0.6 constrains 
x-matter content ilx ^ 0.1 for ~ 0.1, which can fit the observational data by Peacock-Dodds well. Alternatively, if we 
consider the universe to be x-matter dominated, then its sound speed is found to be constrained as < 10~^ at least at the 
epoch of galaxy formation. Although x-matter models are consistent with all current observations, they leave imprints in the 
CMB anisotropy and matter power spectrum that should be detectable in future observations. 
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APPENDIX A: EFFECTIVE POTENTIAL FOR X-MATTER MODEL 

We can construct the effective potential of a scalar field for tlie x-matter model. Tlie metliod of tlie construction is similar to 



that given in Nakamura & Chiba( 

We consider a scalar field (p minimally coupled to gravity. Then the energy density and the pressure of x-matter are given 
in terms of the scalar field as 

Px = (Al) 

Px = i^' + l^W, (A2) 

where V{(j)) is the potential of the scalar field. By the use of Eq.(^ and Eq.(|l|), the scalar field and the potential can be 
written in terms of the scale factor 



r,T// \ 2(c^-Wo) , (l-ci)(l+Wo) -3(l+c2) ..„s 

2V(a) = p^-p^= \. 2 Pxo + ' 2 -P^oa. ^ -\ (A3) 



da J 

is given by the integral 



= (aH^j =(l+wo)Pxoa-'^'-""=^ (A4) 



0(a)- 00 = ± r J 3(l + ci)(l + ^o) du^ 

Then V is given implicitly as a function of (p. The shapes of V{(j>) for — 0.0, 0.1 are shown in Fig. 14 from a = 1 to a = 10~^. 
The scale is arbitrary. Since V{(j)) involves a tiny (but positive finite) constant term, the universe will be dominated by this 
constant term and behave like the cosmological constant dominated universe in future. 



APPENDIX B: FORMULAE FOR A SINGULAR ISOTHERMAL SPHERE LENS 

Here we give useful formulae for a singular isothermal sphere lens which are necessary in the text. 
The deflection angle a by an isothermal sphere lens is given by 

4GM(< b) 
b • 

Here b is the impact parameter and M(< b) is the mass projected within the impact parameter b 



(Bl) 



M(< b)=2Tv I E{r)rdr. (B2) 
Jo 

S(r) is the surface density and is written for a singular isothermal sphere in terms of the velocity dispersion v as 

^ir) = (B3) 
Therefore the deflection angle is written as 
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(B4) 

The lens equation is then 

4> = e~AW^, (B5) 

L>OS 

where cf> and 9 are the angular position of the source and the image, respectively. The critical impact parameter is thus 

, . 2DolDls ,T-,a\ 
bcr = inv — . (B6) 

Dos 



Therefore, the cross section is given by 

,2 ic 3 4 / DqlDls 
V Dos 

If the alignment is close enough such that (j) < q ^^^ , lensing produces multiple images at 



(B7) 



e± = 0±a^, (B8) 
Dos 

An amplification due to lensing is given by 

ede 

A = 



(B9) 

The two images are amplified by factors 

A± = ^ ± 1. (BIO) 

<pDos 

This gives a total amplification 

A = 2^. (BID 
(pDos 

This means that the probability distribution of total amplification for multiply imaged systems is 

P(A)dA = 8A~^dA, for A>2. (B12) 

This paper has been produced using the Royal Astronomical Society/Blackwell Science I^T^^X style file. 



FIGURE CAPTIONS 

Figure 1: The contour plots for Hoto in the {wo,flx) plane. The contour levels are Hoto — 1.0, 0.9, 0.8, 0.7, 0.6(only for 
Cg = 0.15) from top-left to right. 

Figure 2: The lens optical depth for a source at redshift Zs in a flat universe relative to that in the Einstein-de Sitter 
model(f2Af = 1) for {a.)flM = 0.1, (b) Qm = 0.3, (c) Q,m = 0.5, and (d) Q,m = 0.7. Solid curves correspond to — 0, 
dotted ones = 0.1. Each type of curves from top to bottom corresponds to too = —1, —0.8, —0.6, —0.4. Note that wo — — 1 
corresponds to the cosmological constant model irrespective of c^. 

Figure 3: The expected number of lenses for = 0(a),0.05(b),0.1(c),0.15(d) and luo = -1, -0.8, -0.6, -0.4, -0.2,0 as a 
function of ilx. Each curve from top to bottom corresponds to wo = —1, —0.8, —0.6, —0.4, —0.2,0. lOo = —1 corresponds to 



the cosmological constant model irrespective of c^. Data is taken from HST snapshot survey (Maoz et al. 1993) 



Figure 4: Likelihood contours for flat cosmological models projected onto constant Cs plane. Likelihood levels are, from left 
top to right down, 68.3%, 95%, 99%. 

Figure 5: Time evolution of x-matter density fluctuations Ax at the wave number k = O.lMpc for different values of wq and 
C3. Solid and dotted lines are for wq — —0.8, = 0.2 and wo = —0.8, = 0, respectively. A dashed line corresponds to the 
standard CDM model, i.e., wo — and — 0. We take 0,^ ~ 0.95 and Qb ~ 0.05 as cosmological model parameters. 

Figure 6: Jeans scale as a function of redshift. Solid lines correspond to wq = —0.2, while dashed lines to wo ~ —0.8. 



Figure 7: Total matter transfer functions. We take 
© 1998 RAS, MNRAS 000, §-0 



— 0.95 and Qb ~ 0.05 as cosmological model parameters, (a): 



10 T.Chiba et al. 



changing wo with — 0. From left to right, wo = —0.8, —0.6, —0.4 and —0.2. Dashed and dotted lines are CDM models with 
the shape parameter F = 0.5 and 0.25, respectively, (b): changing with wg — 0. From left to right, = 0.15, 0.05, 10"'* and 
10"^ 



Figure 8: CMB angular power spectra. Models are same as Figure 5. 



Figure 9: Contours of as in the wq — plane for various values of Cs- COBE 4yr normalization is employed. We take 
h = 0.5 and Qb ~ 0.05. From top to bottom in each panels, as ~ 0.2, 0.4, 0.6, 0.8 and 1.0 (bold). 



Figure 10: Same as Figure 9 but h = 0.7 and Sis ~ 0.04. 



Figure 11: Total matter transfer functions for models with as = 0.6. Panels (a)-(c) are models with h = 0.5 and Qb = 0.05 
and Panels (d)-(f) are with h = 0.7 and Sis = 0.05. Solid lines are indicated in each panels. Dashed and dotted lines are CDM 
models with the shape parameter F = 0.5 and 0.25, respectively. Since the transfer functions of the model with wo = and 
S7x ~ 0.025 and the one with wo = —0.5 and Six ~ 0.04 are identical in panel (c), we do not plot the later one. 



Figure 12: Matter power spectra for 'best fitted' models. From Figure 11, we chose models with transfer functions similar 
to the ones of F = 0.25. (a): models with h = 0.5 and Sis = 0.05. Solid and dashed lines are models with wq = —0.8 and 
Six = 0.09 and wo = and Six = 0.025. (b): models with h — 0.7 and Sis = 0.04. Solid, dotted and dashed lines are models 
with Wo — —0.8 and Six ~ 0.12, wo — —0.5, Six ~ 0.06 and wo = and Six = 0.03, respectively. The symbols denote the 
observational data by Peacock and Dodds(1994). We multiply 0.8 to these observational values since we are only interested 
in the shape of the matter power spectrum and Peacock and Dodds' data provide roughly 10 — 20% larger as compared to 
the value quoted in the text. 



Figure 13: CMB angular power spectra for models with as ~ 0.6. Models are same as Figure 11. Dotted lines are C^'s of 
CDM models. Namely CDM models with S^o = 1, ^ = 0.5 and Sis ~ 0.05 for panels (a)-(c), and with Slo = 1, /i = 0.7 and 
Sis = 0.04 for panels (d)-(f). 

Figure 14: The effective scalar field potential for x-matter model with — 0,0.1 and wo — —0.8. Scales are arbitrary. 
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